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Abstract. - A continuum Landau theory for the micro-elasticity of membranes is discussed, 
which incorporates a coupling between the bilayer thickness variation and the difference in the 
two monolayers' tilts. This coupling stabilizes a new phase with a rippled micro-structure. 
Interactions among membrane inclusions combine a dilation-induced attraction and a tilt- 
difference-induced repulsion that yield 2D crystal phases, with possible coexistence of different 
lattice spacings for large couplings. Inclusions favoring crystals are those with either a long- 
convex or a short-concave hydrophobic core. 



Lipid molecules in water spontaneously form fluid bilaycrs in which hydrocarbon tails are 
shielded from contact with water JlJ, B. In nature, lipid membranes constitute the walls 
surrounding living cells, and generally host a large number of (protein) inclusions ||. Self- 
assembled surfactant membranes can form various phases, e.g., lamellar (L a ), vesicular (L4), 
or sponge (L3) phases. At low temperatures, the molecules tilt relative to the membrane 
normal, forming the Lpi phase which has an internal degree of freedom similar to that of 
the liquid crystalline smectic-C phase. The coupling between tilt and membrane curvature 
can produce a shape instability yielding the Pp> or "ripple" phase 0), [||- Recently, a new 
degree of freedom has been introduced by Seifert et al.\ a tilt- difference between the two 
membrane monolayers. The coupling between tilt-difference and membrane curvature can 
produce instabilities yielding rippled phases, bilayer tubules and bicontinuous phases H, very 
much like in the case of membranes containing nematogens or anisotropic inclusions Q. 

In this letter, we study the effects of the coupling between tilt-difference and membrane 
dilation in ordinary L a membranes. We find a new type of "ripple" phase that has not 
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Fig. 1. - Lipid membrane subject to the ripple instability produced by the coupling c m ■ Vu between 
the dilation (u) and tilt-difference (m) modes. 



yet been evidenced. The tilt-difference and dilation modes are generically excited by protein 
inclusions having a convex or a concave shape, or a hydrophobic core with a thickness different 
from that of the bilayer. We calculate and discuss the short-range interactions between such 
hosts, which form the most general up-down symmetric inclusions. 

At lengths scales large compared with molecular dimensions, membranes are traditionally 
described by the shape h(r) of their midsurface, with r = (x,y). Their elasticity is governed 
by the Helfrich curvature Hamiltonian |(k + k)(V 2 /i) 2 — |k (didjh)(didjh) Q. There is no 
(V/i) 2 surface tension term: membranes are self-assembled systems that optimize their area 
per molecule (or their coarse-grained area at larger scales.) We use here a different approach 
since we are interested in the microscopic elasticity of the bilayer pl]|-Pq|. To simplify, we only 
consider membranes perfectly symmetric with respect to their midsurface (a general theory 
treating the two monolayers independently will be presented elsewhere [|l0|.) We develop a 
phcnomcnological Landau theory for the membrane dilation u(r) (thickness variation), and 
the molecular tilt-difference m(r). The latter is defined as half the sum of the projections 
onto the (x, y) plane of the unit vectors, oriented tail-to-head, parallel to the molecules in the 
two monolayers (fig. |l|). The model free energy involves all quadratic terms and first order 
derivatives that satisfy rotational symmetry: 

f=-Bu 2 + ^A(Vm) 2 + cm- Vu+ -t' rh 2 + -K[ (V • m) 2 + -K' 2 (V x m) 2 . (1) 
2 2 2 2 2 

Due to the tendency of the molecules to orient perpendicular to the chain-water interface, 
we expect c > 0. For biological membranes, the typical energy and length scales are given 
by k — 25fc B T (KT 12 erg) and £ - 20 A(monolayer thickness) 0. This yields B « k/£$ ~ 
6 x 10 14 erg cm which agrees with the experimental value of the area-stretching coefficient 
B(2£o) 2 ~ 100 erg cm" 2 jjj. The dilation term \Bu 2 has contributions from the water- 
chain (oil) surface tension and from the polymer-like stretching of the chains. The dilation- 
gradient term iA(Vu) 2 therefore originates from the extra cost of modulating the stretching 
of the chains and increasing the water-chain area. Hence, we expect A ~ B^ 2 ~ 25 erg cm" 2 , 
dimcnsionally, whether or not the (macroscopic, effective) tension of the membrane vanishes. 



On this point, our model disagrees with those of Dan et al. 18 . On the same basis, the other 
coefficients c, t' are expected to be ~ k/£ 2 and the K'^s are expected to be ss k. Neglecting 
higher-order gradient terms in (^) is in fact only justified when B is small, which is the case 
in the vicinity of a transition to a more ordered Lp or Lpi phase. For a microscopic y/X/B, 
we expect nonetheless to get qualitatively correct results. 
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Fig. 2. - Protein inclusions with thickness mismatch Uo and tilt-difference angle ao- Inclusions forming 
2D crystals have either a long-convex (a) or a short-concave (b) hydrophobic core, (c) Wigner-Seitz 
cell approximated by a circle in an array of inclusions. 



"Ripple " phase. - The linear stability of the hat membrane can be studied in the reciprocal 
space, where the free-energy density per mode associated with (w) can be diagonalized to 
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where 



and 



are the projections of the tilt-difference Fourier component rhk parallel 



and perpendicular to k, respectively, and B(k) = B + \k 2 > 0, t'(k) = t' + K' 2 k > 0, 
JV(fc) = (f + K[k 2 ){B + \k 2 ) - c 2 k 2 . The stability of the flat membrane is dictated by the 
sign of N(k). When the latter is positive, the minimum energy corresponds to uu = rhk = 
and the flat membrane is stable. Defining the dilation-, the tilt-difference- and the coupling- 
characteristic lengths, as 



(3) 



1 < 0. Therefore the 




respectively, N(k) < is equivalent to £ 2 £' 2 fc 4 + (£ 2 
instability occurs for £ > (£ + £')/2 at a nonzero wavevector k c — (££') 
tilt-difference modulations being in quadrature (fig. [j]). 



-1/2 



the dilation and 



Short-range interactions among inclusions. - Understanding the membrane-mediated inter- 
actions among inclusions has recently attracted much interest. Conical inclusions coupling 
to the local membrane curvature are subject to long-range elastic and Casimir (fluctuation) 
forces |22| that add up to the standard direct forces. On the other hand, short-range inter- 
actions arise from local structural perturbations: integral proteins have a central hydrophobic 
region spanning the hydrophobic core of the membrane and two polar extremities protruding 
outside. Any thickness mismatch between the protein and the bilayer hydrophobic regions 
results in a thickness perturbation that yields membrane- mediated interactions [p~2| -| 16 . 

Here, we focus on the effects of the tilt-difference distortions that are naturally excited 
by inclusions with a convex or concave hydrophobic region. To simplify, we assume £' = £ 
and hence £ < £ for the membrane's stability, and we restrict attention to inclusions having 
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Fig. 3. - (a) Sketch of the membrane structure around an isolated inclusion with boundary dilation 
(u) and tilt-difference (a) normalized to 1 (ro/£ = 3, x = —2, <f> — 0.75 x 7r/2). (b) Interaction energy 
F in units of 10 ksT vs. separation for zero dilation-tilt-difference coupling (ro/£ = 3, x = 0.5, = 0). 



revolution and up-down symmetry. We denote by uo their hydrophobic thickness mismatch 
and by ao their tilt-difference angle (fig. |]a,b). Their mean field interaction can be derived from 
the energy of the membrane equilibrium structure, which obeys the Euler-Lagrange equations 
associated with ([!]) , 

Bii-AV 2 « = cV-m, (4) 
]fiV(V'm)+^Vx(Vxm) = -cVu. (5) 



f'm- 



To study the physics of a collection of inclusions we consider the Wigner-Seitz cell surrounding 
one inclusion with radius ro in a lattice (fig. ^p) and we approximate it by a circle of radius 
R jt2| , jt4| . The problem therefore acquires revolution symmetry. In an hexagonal lattice, this 
amounts to neglecting high-order Fourier harmonics (6, 12, etc.), and in a gas to considering 
that the first neighbors effectively screen the other inclusions. Assuming radial symmetry, 
u = u(r) and m = a(r) f , the equilibrium solutions are given by the real part of 



u{r) 
a(r) 
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(6) 
(7) 



where i = y— T, the I's and the K's are modified Bessel function, sin0 = (1 < £), and Ai, 
A2 are two dimensionlcss complex constants that we determine form the boundary conditions: 



u ■ 



a 



u'(r)\ R = 0, a\ R = Q. 



(8) 



The distortion energy F stored in the region ro < r < R p3[ can be transformed by integration 
by parts into F = n [Xruu' + cruet + K' x (ra! + a)a] r , which can be scaled to 



ttB r Q £ ul 



F [x,<t>. 
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(9) 



For typical integral proteins, we expect ro — 3 £ with £ ~ 20 A. Both the dilation and the 
tilt-difference distortion around an inclusions relax exponentially to zero on a distance ~ 100 A 
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Fig. 4. - (a) Interaction energy F in units of 10 k^T vs. separation for nonzero dilation-tilt-difference 
coupling (ro/£ = 3, x = —3, (j> = 0.75 x n/2). (b) Phase diagram for a membrane with £ ~ 20 A 
hosting inclusions with radius ro — 3£. (D) disordered. (K) crystal. (K n ) region where n distinct 
stable or metastable crystalline phases are possible. For x ~ 0, the crystal spacing is so small that 
the inclusions effectively aggregate. 

(~ 5£), as shown in fig. ||a. The coupling between dilation and tilt-difference generates 
oscillations and overshoots. For zero coupling, F can be splitted into a dilation-induced 
attraction F u and a tilt-difference-induced repulsion F a that combine to yield a minimum for 
small enough x (fig. ^p). With the values of the material constants previously discussed and 
typically uq ~ 0.2 £ (~ 4 A), we expect ctQ ~ x x 10 deg, and ttB ro £ u\ ~ 10 k^,T, which sets 
the magnitude of the interaction energy, and in particular the depth of the well in fig. |3|b. For 
large couplings, interferences between the membrane oscillations yield several energy minima 
(see the two wells of depths ~ 25 k^,T and 3 k^T in fig. ^a). 

Figure ^b shows a typical phase diagram for a collection of up-down symmetric inclusions. 
The criterion for a crystal phase is the existence of a well deeper than k^,T. Solid lines indicate 
first-order phase transitions; dashed lines separate regions where different metastable crystals 
with different lattice spacings can exist. Remarkably, crystal phases occur preferentially for 
x < 0, i.e., for inclusions such as depicted in fig. |^a,b. 

Discussion. - The instability produced by the coupling between tilt-difference and dilation 
defines a new "ripple" phase of membranes. The corresponding corrugation actually forms 
a micro- structure, since the undulation period compares with the bilayer thickness (except 
maybe in the vicinity of a Lp or Lp> phase). The amplitude of the corrugation, governed by 
the bilayer dilation elasticity, will probably not exceed a few angstroms: these ripples will be 
difficult to detect by electron microscopy or STM techniques. The present phase cannot be 
mistaken with the Ppi phase, which exhibits much larger (height) amplitudes ~ 45 A p4| . 

Up-down symmetric membrane inclusions generally have a slightly convex or concave hy- 
drophobic core of thickness different from that of the bilayer. Due to the strong effective 
attraction between hydrophobic parts, such inclusions will excite the coupled tilt-difference 
and dilatation modes, which in turn will mediate short-range interactions between them. The 
thickness mismatch creates an energetic dilation corona around the inclusions and yields an 
attraction between like inclusions: no extra distortion occurs when the coronas overlap since 
the boundary dilations match. The tilt-difference, however, yields a repulsion between like 
inclusions: going from ao to — c*q, it develops a strong gradient when the coronas overlap. 
Inclusion producing no tilt-difference aggregate, while inclusions producing a nonzero tilt- 
difference either repel one another or form 2D crystals. The latter situation arises for small 
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tilt-differences, or when the dilation corona extends further than the tilt-difference corona 
(£ > £'). In both cases, the attraction which dominates at "large" distances is overcome by 
the repulsion at short distances. 

When the coupling between dilation and tilt-difference is large, the distortions in the coronas 
exhibit damped oscillations and the intcrparticle potential develops several minima. This 
implies the possible coexistence of different crystals of inclusions having different lattice spac- 
ings, and also the possibility of low density crystals (separation between inclusions' boundaries 
~ 180 A for the secondary minimum in fig. ^). The inclusions most likely to form 2D crystals 
are those with either a long-convex or a short-concave hydrophobic core, i.e., those disfavored 
from the point of view of the m • Vu coupling. This is because the gradient of u being 
more costly, the dilation corona extends (favoring "long-range" attraction), while due to the 
effective shift in the tilt-difference's quadratic potential, the dilation corona shrinks (making 
the repulsion occur only for smaller separations). Conversely, short-convex and long-concave 
inclusions have a dominant repulsion and should form disordered phases. 

In replacing the Wigner-Seitz cell by a circle, we have neglected high-order Fourier harmon- 
ics. Taking them into account would probably change the detail of the interaction potential, 
but not its essential features. Real inclusions are not up-down symmetric: their average 
conical shape effectively adds a monotonic repulsion [ fl0| , which can be neglected provided the 
corresponding angle is small compared to the tilt-difference angle. 
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